Analytical theory of dark nonlocal solitons 
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We investigate properties of dark solitons in nonlocal materials with an arbitrary degree 
of nonlocality. We employ the variational technique and describe the dark solitons, for the 
CN ' first time, in the whole range of degree of nonlocality. © 2010 Optical Society of America 
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n , . Spatial dark optical solitons are localized light intensity dips in the infinite constant background [T}]3]. 

r^ , Like their bright counterparts, they propagate preserving their spatial profile. This is a result of a balance 

between diffraction which tends to broaden their size and nonlinearity of the medium which counteracts this 

process. The existence of dark spatial solitons requires the nonlinearity of the medium to be negative, or 

self-defocusing, i.e. refractive index of the medium must decrease with light intensity [T]. In the simplest case 

^ . of Kerr medium the change of the refractive index is just proportional to the light intensity. In this special 

'/^ , situation and in one-dimensional case the soliton amplitude u{x, z) is described analytically by the following 

U-v \ relation 

0\J ' u{x, z) — B tanh(x — vz) + lA, (1) 

fm^ . where x and z denote the transverse and propagation coordinates, respectively. B and A are constant and 

f^ , V is soliton transverse velocity, v = A/ B. In case of more complex relation between light intensity and the 

nonlinearity the solitons profiles, except for some special models, can be only found numerically [5]. All the 
typical nonlinear models are local i.e. the nonlinear response in a particular point is determined solely by the 
S^ ■ light intensity in the same point. However, there has been recently strong interest in the so called nonlocal 

nonlinearities [4]. In those models the relation between nonlinear response and the intensity is spatially 
nonlocal. Physically it means that the nonlinearity in a particular spatial location is determined by the 
light intensity in a certain neighborhood of this location. Typical nonlocal systems involve, either transport 
processes such as heat [S] or ballistic atomic transport [5] and diffusion [7] , charge separation [5] or long-range 
interaction as in dipolar Bose Einstein condensate [9] or nematic liquid crystals |10| . Nonlocal nonlinearity 
has been shown to have a stabilizing effect on nonlinear structures [4l [TTm3] enabling formation of complex 
solitonic entities such as multi-pole solitons or azimuthons |14j . It also affects soliton interaction inducing 
a long-range attractive forces between solitons [T5H17) . In case of self-defocusing nonlinearities, nonlocality 
has been shown to modify soliton interaction leading to attraction of dark solitons and formation of their 
bound states [T8l[T9]. Unlike the Kerr nonlinear models where the solitons can be described analytically, 
fully nonlocal models with arbitrary degree of nonlocality can only be treated numerically [20 [ [21 ] . The only 
exceptions arc the special case of the so called weak and strong nonlocality, when the steady state bright 
and dark soliton could be found in analytical form j^l^ . 
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In this work we describe analytically, for the first time, properties of dark nonlocal solitons. To this end we 
employ a variational approach and show that it enables one to retrieve the major features of dark solitons 
in a general nonlocal regime. 

The evolution of one-dimensional optical beam with an amplitude u{x, z) in nonlocal defocusing medium 
is governed by the following nonlocal nonlinear Schrodinger equation (NLS) [5] 



du 1 d'^u f^'^ 

*^+2a^-"y_^ ^(^-C)K,^)l'rfe = 0. (2) 

Here we used the phcnomcnological model of the nonlocal nonlinearity with R{x) being the nonlocal response 
function. Its width determines the degree of nonlocality [35]. In particular, for R{x) ~ S{x) the above 
equation describes standard Kerr local medium. While this is only a phenomenological model it nevertheless 
describes very well general properties of the nonlocal media. Moreover, for certain form of the nonlocal 
response function this model represents the actual physical system. For instance, this is exactly the case of 
nematic liquid crystals, when the long-range interaction-mediated nonlinear response function, under certain 
approximations [23], can be written as R{x) oc exp (— |x|/cr) with a being the degree of nonlocality. Also, the 
same exponential response describes nonlinear interaction in quadratic media [24j . 

To analyze the nonlocal NLS equation, we will employ the Lagrangian approach which in case of dark 
solitons has been formulated in |25| . Following this work wc find that the renormalized Lagrangian density 
corresponding to the NLS Eq.Q is given in the following form 
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-l(iup-i) / i?(x-e)(K^,^)r-iK. (3) 



To proceed further we need to specify the nonlocality. For the sake of simplicity and analytical tractabil- 
ity and without the loss of generality we consider here the rectangular profile for the nonlocal response 
function [22] R{x): 

R{x) = [ ^ -^-^-' , (4) 

I otherwise. 

with a defining the degree of nonlocality. For cr — > we get R{x) — > S{x) and the model describes just the 
local Kerr nonlinearity. The accuracy of the Lagrangian approach depends on the functional choice of the 
variational solution. Here we use the obvious ansatz 



u{x,z) ~ Bt'dn]i[D{x — xq)] + iA, 
A'^+B^ = 1, 



(5) 



which represents an exact soliton solution in the local regime. All parameters A,B,D, and xq are assumed 
to be functions of the propagation variable z. Substituting Eq.(l5]) and Eq.Q into the Lagrangian density 
Eq.(j31) and integrating over x, we get the "averaged" Lagrangian, 
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C{u)dx 



, dxo 
dz 



'AB + t-Ai\ ^ 



- ^B^D 



B^ 
D 



csc\?{Da) - 



Then wc find the corresponding Euler-Lagrangc equations as 



coth(L» cr) 



(6) 



coth{Da) 



Da 



D^ 



(7^csch^{Da) 



dB _ 
dz 

1 
" 3S2' 



dxo 
dz 



= A 



D B f n, , coth Do- 

csch^ Do- 5^ - 

3B D \ ^ ^ Da 



(7) 



Notice that for a ~ the above formulas give B = D = const and dxo/dz ~ A which are the exact sohton 
parameters for local Kerr solitons [25] . The formulas Eq.© constitute the main result of this paper. They 
give, for the first time, the analytical relation among parameters of the dark nonlocal soliton, for arbitrary 
degree of nonlocality. In particular, one can use it to show the dependence of the soliton width (oc D~^) 
as a function of the degree of nonlocality (a). This relation is depicted in Fig.l by the solid line. Here the 
soliton width is normalized to its value in the local regime. The graph clearly shows that the width of dark 
soliton is a nonmonotonic function of the nonlocality. It decreases first for small ct, reaches its minimum and 
then monotonically increases. This nontrivial analytical result agrees with that found earlier in numerical 
simulations [see Fig. 3 in [H]]. On the physical grounds the initial decrease of the soliton width can be 
explained by the fact that weak nonlocality causes the nonlinear index change to advance towards regions 
of lower light intensity. As a result the waveguide induced by the soliton becomes slightly narrower and so 
does the soliton. Situation becomes different for high degree of nonlocality. For large a the refractive index 
modulation expressed by the convolution integral in Eq.Q becomes weaker and broader, acquiring wide 
rectangular profile and resulting in increased width of the soliton. 

It is instructive to compare the above variational calculations with the exact analytical solutions which 
can be obtained in two limiting regimes of weak and strong nonlocality. In the former case when the response 
function is much narrower than the soliton width, the convolution term in Eq. ^ can be expanded in Taylor 
series leading to the following form of the nonlinear response 

R{x - 0\ui^, z)\'d^ « \u{x)\' + 7 '^y , (8) 

where we introduced parameter 7 = <t^/6. The NLS with such nonlinear term can be solved analytically as 
shown in Rcf. 122] . The weakly nonlocal limit can be recovered from the general variational solutions Eq.([7|) 
by expanding it into a Taylor series around ct = 0. After retaining only the most significant terms we obtain 

^ = 0, 
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which coincides exactly with the equations derived directly from the Lagrangian representation of the weakly 
nonlocal nonlinear Schrodinger equation j26] The exactly found width of dark soliton in a weakly nonlocal 
regime solution is shown in Fig.l by dashed (red) line. It is evident that the agreement between variational 
and exact solutions is indeed very good. This is an obvious consequence of the fact that the ansatz, Eq.([5|), 
is actually pretty close to the exact solution in this regime. Different situation arises in the so called highly 
nonlocal limit. It has been shown that in such limit the nonlocal NLS becomes linear with the convolution 
term being just proportional to the response function. In this regime this equation can be solved exactly. 
For the rectangular response function the solution is 

u{x) ^ B sin ('^^^^^J^\ + iA (10) 

This gives monotonic increase of the soliton width with nonlocality D^^ ex a. This of course differs from the 
variational solution Eq.© which for large degree of nonlocality gives D^^ oc a^'^ . This discrepancy is most 
likely caused by the inadequacy of our ansatz in the highly nonlocal regime as will be evident below. 

In Fig. 2 we illustrate the effect of nonlocality on the formation and propagation of single dark spatial 
soliton. We used split step Fourier method to integrate numerically the nonlinear Schrodinger equation. Fig. 2 
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Fig. 1. Width of the dark nonlocal soliton (normalized to its value in the local limit) as a 
function of the degree of nonlocality a. solid (blue) line - variational results with rectangular 
response function; dashed (red) line - exact solution in a weakly nonlocal limit Ref. |22) . 
The inset depicts the magnified region of weak nonlocality. Here i? « 1, ^ w 
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Fig. 2. Dynamics of formation of dark spatial solitons in (a) almost local, a = 0.2 and (b) 
nonlocal, cr = 3, media with rectangular response function. The insets depict the soliton 
profile at the beginning (bottom) and end (top) of the propagation distance. 



depict dynamics of dark soliton excited by the initial conditions represented by the variationally calculated 
profiles with ct = 3. Graphs (a) and (b) correspond to almost local (cr = 0.2) and nonlocal (cr = 3) cases, 
respectively. It is clear that nonlocality tends to expand the width of soliton. The strong dispersive waves 
visible in graph (b) reflect the fact that the initial conditions do not correspond to the stationary soliton 
profile. This is obvious after noticing that unlike the local case when the soliton intensity monotonically 
grows from its minimum to its background value, the nonlocal intensity profile almost always exhibit strong 
overshot [clearly visible in the inset in Fig. 2(b)]. Unfortunately, the ansatz Eq.(l5]) does not reflect this 
property and hence excitation of dark soliton using our variational solutions is always accompanied by the 
emission of diffractive waves. 

In conclusion, wc studied properties of single dark solitons in nonlocal nonlinear media. We used the 
variational approach an the rectangular model of the nonlocality to derive, for the first time, analytical 
relations for soliton parameters for arbitrary degree of nonlocality. We showed that this approach, while 
approximate, faithfully represent properties of dark nonlocal solitons. 
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